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ABSTRACT 



We give one more proof in two and three space dimensions that the irregular solution of the 
Schrodinger equation, for zero angular momentum, is in fact the solution of an equation con- 
taining an extra "delta function". We propose another criterium to eliminate the irregular 
solution which is to require the validity of the virial theorem of which we give a general proof 
in the classical and quantum cases. 



1 Introduction 

These pedagogical remarks were stimulated by an exchange of ideas with Yves Cantelaube. 

We re-examine the following question: in the case of a central potential why choose the 
"regular solution" at the origin of the Schrodinger equation, while the irregular solution may 
also be square integrable? It is indeed easy to show that for any negative energy there exists, 
under weak conditions on the potential, a solution exponentially decreasing at infinity which, 
outside the conventional eigenvalues, is irregular at the origin and square integrable. This is 
explicitly demonstrated in Section 2. I would just like first to present, in my own way, the 
standard argument (Diracfl], Van Hove [2], Cohen- Tannoudji et al.[3], Basdevant et al.[4]) that 
the irregular solution is not a solution of the original Schrodinger equation but a solution of an 
equation containing an extra delta function. Then, I would like to propose another criterium: 
the violation of the virial theorem by the irregular solution. I refuse to enter into the questions of 
self-adjointness. You prove that a certain operator is self-adjoint on a certain domain. Choosing 
the domain is more or less choosing the answer. 
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2 Proof of the existence of singular, normalizable solu- 
tions 



Consider the reduced Schrodinger equation in three dimensions for i = 0: 

^ + V(r)-Ey = 

u = rip. (1) 
Assume V — > for r — > oo. Then for £ < 0, there exists i? such that 

V(r) - E > K 2 > for r > R. (2) 

Consider the solution such that 

u\R) = 0. 

By comparing Eq. (1) with 
it is easy to see that 

u(r) > u(R) coshK(r - R) 

for r > R, (4) 

and that u(r) is increasing. 

We can define another solution of Eq. (1), v which has a constant Wronskian with u: 

u'v — v'u = 1 (5) 

A particular solution of (5) is 

v = u —-— (6) 

Jr u 2 (r) 



From (4) we get 

dr' _ 2 exp - K(r - R) 

lr u(R)coshK(r - R) ~ K u(R) 
The solution v can be continued for r < R till the origin. 



f 00 dr' _ 2 exp — K [r — K) 

V< Jr u(R)coshK(r - R) ~ K M) ^ ' 



Equation (1) can be written in integral form as 

rRo 



v = v(R ) + (r - Ro)v'(Ro) - (r' - r)(V(r') - E)v{r')dr' (8) 

J r 

Hence, for r < R , 



M{r,R ) < \v(Ro)\+Ro\v'(Ro) 
fRo 

+ 

/o 



rRo 

/ r'[\V(r')\ + \E\]M{r,R )dr', 
Jo 



where 

M(r,R ) = sup\v(r')\. 

r <r' < Rq 

If 

/ r'\V{r')\dr' < oo, (9) 
Jo 

we can choose R in such a way that 



/ r'[\V(r')\ + \E\}dr' <e<l 
Jo 



and hence 



MW )<M^W (10) 



so |t> (r)| has a fixed bound, independent of r. 

Condition (9) is only sufficient. In particular one can have acceptable oscillating potentials 
as shown by Chadan and collaborators [5] which violate (9). 

Now, since the wave function is given by 

ip = const, rv , 

it is clear that ip is normalizable because / \vfj\ 2 d 3 r converges both at the origin and at infinity. 
This is valid for any E < 0, in particular for a value E which is not a conventional eigenvalue 
(for which v would tend to zero at the origin). 



3 The delta function in the radial case with three-dimensions 

We want to prove that the solution of the Schrodinger equation 

4> = -u (11) 
r 

cannot behave like £ at the origin. 
Take 

/ = -r^== 2 (12) 

A/ = — — r 2 — / = -3 ^ = -g (13) 
r 2 dr dr (r 2 + a 2 ) 5 / 2 



Let us prove that for a — > 0, g tends to 4nS 3 / (r). 
1) for r = 0, g — > oo for a — > 
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2) for any r > 0, g — > 0, for a — > 

3) /^(rl=4 7 r/ oo I |g|^ = 47r. 

Hence for 

a -> 0, p -> 47T<5 3 (f) 

A 3 (--) = +4vr<5 3 (f). (14) 
r 

Notice that the sign is obvious since — 1/r is a subharmonic function. Hence, as announced 
the irregular solution of the Schrodinger equation is not acceptable. 



4 The case of two-dimensions 

The regular solution for zero azimuthal angular momentum behaves like 1, the irregular solution 
behaves like lnr. We use 

/ = In Vr 2 + a 2 = - ln(r 2 + a 2 ) (15) 

Af = ~4-r^- \n(r 2 + a 2 ) 
J 2rdr dr y ' 

2a 2 

9 (16) 



(r 2 + a 2 ) 2 

Let us prove that g is proportional to a delta function in the limit a — > 

1) again g — *ooforr = 0a^0 

2) -> for r > a -> 
3) 

/ rf 2 ™ — O.ir I rdr 

(r 2 + a 2 ) 
Hence 

^ - 2m5 2 (f) (18) 

for a — > 0. Again the sign is correct because lnr is subharmonic since lnr tends to MINUS 
infinity for r — > 0. 



//■ 2(2 2 
d ^ = 2 "J rdr (^^ = 2 " (17 > 
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5 The virial theorem 



The virial theorem should hold classically and "quantically" . Let me remind the classical proof 
in the general case. 

If Fi = — VjV(xi, x 2 , x n ) , (19) 

Newton's equations are 

m,^ = -V,y . (20) 

E^-^ = -E^w (2i) 



Hence 



and 



" 1 S ^ • = " S m ^ |o + / E ^ ("f J rft = / E aV 4 ™ (22) 
So if you deal with a confined system (big problem!) 

lim — f y^mjfj^ | = . 
T^oo dt J 

Strictly speaking, we know only that the above parenthesis is bounded for almost every T if 
the X{ S £1X6 bounded. This inconvenient could be avoided by smoothing the cut-off at T in the 
integrals. 

With this caveat, the average over time of the kinetic energy is: 

If some particles escape at infinity with velocities Vi (nobody knows if this will happen in the 
solar system), we get 

li m I [ T IVm l f^V^ = E-^+ lim - f dt (24) 

t^oo T Jo 2^ \dt j ^ 2 1 1 t^oo T Jo 2 v ; 

and hence we have an inequality instead of an equality. 

In the quantum case, the average over time of the kinetic energy is replaced by an average 
over space of the kinetic energy operator. For the radial case the virial theorem can be proved by 
hand, but, following Thirring[6], the best general proof uses scaling. Consider the Hamiltonian 

H = -'£^ + V(x 1 ,x 2 ,...,x n ) , H^ k = E k ^ k (25) 



if we change the scale x« — > Xxi, the energy will be unchanged for 

A 2 ^ 2m. 



#a = ~ E P~ + ^i, •••> Xx n) (26) 
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so if we consider the Hamiltonian 

\ 2 H x = -'£ 1 ^- + \ 2 V(\x 1 ,...,\x N ) , (27) 

the energy levels e k are 

e fc (A) = \ 2 E k (28) 
so, for A = 1, by the Feynman-Hellman theorem, 

^± = 2E k = 2(V(x 1} x N )) + (mViV(xi, .., x N )) , (29) 

so 

E k = (V(x 1 ,...,x N )) + ^(mV i V) 
and since E k = (T) + (V), T being the kinetic energy operator, 

(T> = ^(Ef i V i >. (30) 

Remark: the choice of the origin is irrelevant because (a ■ VV) = if a is a fixed vector. This 
can be verified in an elementary way, using the Schrodinger equation. 



6 Virial theorem in the three-dimensional quantum case 



First, we would like to give also the "pedestrian" proof of the virial theorem in the 3 dimensions 
radial case, to show the importance of the prescription u — 0. 



Take the Schrodinger equation 

-u n +Vu = Eu , 
multiply by u and integrate. This gives 



-uu 



-J u' 2 dr + J Vu 2 dr = E J 



u dr 



multiply also by 2ru' and integrate. This gives 



— ru 



+ 



r 

Jo 



u' 2 dr + rVu' 



-F 

Jo 



\r 2 dV 2 \ J 

Vu + r——u dr 

dr 



— Eru' 



- E \ u 2 dr 



Assume that u and u' vanish at r — > oo. Adding up the 2 equations, we get 

roo dV 



lim (uu' + ru' 2 - rVu 2 ^j + 2 j u' 2 dr = J 



r— — u dr 

dr 



(31) 



(32) 



(33) 



(34) 



If u(r) ~ r for r — > and if -u'(O) is finite, we get the virial theorem, if w(0) = 1 and if u'(0) 
is not zero, the virial theorem is violated. If V(r')dr' — > ±00 for r — > 0, it is impossible to 
have 1/(0) = because |w'(V)| — > 00 for r — > 0, and this dominates limrV. 

If V is integrable near the origin, one can save the virial theorem by taking u'(0) = 0. This, 
however, is equivalent to looking at a symmetric one dimensional potential V(x) = V(— x) and 
looking at the even levels. Notice that any other boundary condition at the origin violates the 
virial theorem. 

Now the point is that if the virial theorem holds at the classical level, it should hold too at 
the quantum level, according to the sacro-sanct "correspondence principle" of Bohr. 

In the case of Coulomb interactions, we see that 

Id 11 

-r—V[r) = 

2 dr y ' 2 r 

so 

(T) = ~{V) = -E. 

Taking a solution irregular at the origin for the case of Coulomb in 3 dimensions gives 

r 1 1 

(V) ~ - / r 2 dr- x — = -00 
J r r l 

and the virial theorem is completely violated. The same is true for potentials more singular 
than Coulomb. 

For the radial case, we have seen that we can save formally the virial theorem if V is 
integrable at the origin and u(0) = 1 only if u'(0) = 0. This corresponds to an unphysical 
self-adjoint extension. Indeed while / u' 2 dr is finite, / \ Vtp\ 2 d 2 r diverges. 



7 Virial theorem in the two-dimensional quantum case 

In the radial case, the possible behaviour of the wave function at the origin are 

ip ~ c , ip ~ lnr (35) 

If — > c, the kinetic energy is given by 



(T) 



dr 2 = - J i)^i\)d 2 r = ^J (r^p\ ^ 2 d 2 r , (36) 



and the virial theorem is satisfied. More specifically, using the reduced wave function 



u = V2^ , (37) 
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which satisfies 



-u" - ^- + Vu = Eu 



and using the same stragegy as in the three-dimensional case, we get 



(T) 



- lim 

2 r^o r 



\u{r)' 



+ 



rdr = - 
Ar 2 J 2 



dV 
r — — u dr 

dr 



If on the other hand, ip ~ In r near the origin, Vip ~ £ and hence 

2 



/ 



rdr 



It- 



oo 



(38) 



(39) 



(40) 



Therefore, in all cases , for any potential, there is a violation of the virial theorem by the 
irregular solution. This, we believe, settles controversies on the choice of boundary conditions 
at the origin in a paper on the number of bound states in one and two space dimensions [7]. 
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